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Abstract 

The problem of extracting information on new and known N* resonances by fitting 
isobar models to photonuclear data is addressed. A new fitting strategy, incorpo- 
rating a genetic algorithm, is outlined. As an example, the method is applied to 
a typical tree-level analysis of published p{^,K + ) A data. It is shown that, within 
the limitations of this tree-level analysis, a resonance in addition to the known set 
is required to obtain a reasonable fit. An additional P\\ resonance, with a mass of 
about 1.9 GeV, gives the best agreement with the published data, but additional S\\ 
or D13 resonances cannot be ruled out. Our genetic algorithm method predicts that 
photon beam asymmetry and double polarization p (7, K + ) A measurements should 
provide the most sensitive information with respect to missing resonances. 
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1 Introduction 



The spectroscopy of baryons continues to be a subject of great interest in 
intermediate energy nuclear physics, since it an essential component in under- 
pinning our knowledge of the substructure of the nucleon. Central to this is 
the issue of which resonances exist. Constituent quark models such as those 
proposed by Capstick and Roberts [1] predict a large number of resonances 
which thus far have not been observed, whilst for models which limit quark 
degrees of freedom [2], the number of "missing" resonances is less. 

The majority of resonance information has been gleaned from analyses of single 
pion production reactions [3], and the suggestion that the missing resonances 
may couple more strongly to other channels [1,4] has led to a number of ex- 
periments being proposed and carried out at intermediate energy accelerator 
facilities. Of particular interest is the possibility of studying strangeness pro- 
duction, since this opens up the possibility of extracting extra N* information. 

The use of effective field theories is necessary in the resonance region since 
QCD cannot be solved perturbatively at this energy scale. Constituent quark 
models are able to predict quantities such as coupling constants and electro- 
magnetic form factors, but it is not straightforward for these models to predict 
real observables, as a fundamental understanding of the underlying dynamics 
is missing. By using an effective field theory to extract coupling constants from 
fits to data, one can obtain numbers to compare with quark model predictions. 

A complete understanding of the physics underlying the photonucleon data 
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in the resonance region will in principle only be possible in a framework 
that contains all participating processes which include coupling to resonances 
[5,6,7,8,9]. This means taking into account meson production reactions such 
as 7iV — > nN, ttttN, rjN, uN, KA, as well as Compton scattering and 

meson-induced production reactions. This is clearly an enormous task. 

In a model adopting effective degrees-of-freedom, the field associated with each 
resonance has a number of free parameters which are usually determined by 
fitting model calculations to data. In a full-blown coupled-channels approach, 
the number of free parameters could easily end up being over 100. For such 
a procedure to be possible, a sizeable data set for each of the contributing 
channels is required. One difficulty, which is often overlooked, is the process 
of how to obtain a set of parameters which leads to the best description of the 
data. Related to this is the issue of how to assign error bars and confidence 
levels to the values of the extracted resonance parameters. 

An analysis of a single channel, whilst not complete, offers the possibility of 
simplifying the problem by reducing the number of free parameters to a man- 
ageable size, and in being able to identify the most important features. Recent 
attempts to analyse the p(j, K + )A channel in such a framework have high- 
lighted the tantalising prospect that previously undiscovered resonances may 
reveal themselves in the mechanisms involving strangeness production. Mart 
and Bennhold [10] used their tree-level hadrodynamical model to show that 
the total cross-section measured at SAPHIR [11] could be reproduced by intro- 
ducing a _Di3(1895) resonance. This resonance does not appear in the Particle 
Data Group baryon summary table as it had not been conclusively observed, 
but was included in [10] because of the predicted coupling [1] to strange chan- 
nels. In addition to this, a recent re-analysis of pion photoproduction data 
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[12] appeared to find evidence for it. On the other hand, an analysis by Saghai 
[13] argued that by tuning the background processes involved, the need for 
the extra resonance was removed. Both approaches were based on the analysis 
of the SAPHIR data set, which at that point amounted to around 100 data 
points. 

In a previous article [14], we highlighted the problem of extracting reliable 
resonance information from the limited set of p(7, K + )A data points by per- 
forming many independent fitting calculations. Using a fitting strategy based 
on a genetic algorithm we illustrated that a large number of different solu- 
tions resulted in similar x 2 values. We were able to group these solutions into 
two main sets: one with _Di 3 coupling constants close to zero (indicating its 
non-existence), and one with significantly non-zero D Vi coupling constants 
(indicating its existence). We further showed that the measurement of polar- 
ization observables, such as photon beam asymmetry, could have a pivotal role 
in removing ambiguities. 

We have now performed a systematic study into the feasibility of determining 
the combination of resonances which contribute to the ^(7, K + )A reaction. To 
this end, we have used a typical tree-level description of the reaction process 
in combination with a minimization procedure which is specifically designed 
to quantify the relative success of different combinations of resonances. This 
procedure employs a genetic algorithm which is able to search efficiently a large 
parameter space, coupled with a conventional minimization routine to ensure 
convergence. We have benefitted from adding new photon beam asymmetry 
measurements from SPring-8 [15], as well as fitting electroproduction data 
measured at Jefferson Lab [16]. 
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The structure of this paper is as follows. We outline the features of a typ- 
ical tree-level p{^,K + )A reaction model in section 2. Section 3 describes a 
framework which enables the extraction of resonance information from exist- 
ing photonucleon data, and discusses how to evaluate the reliability of the 
information. As an example, the methodology is applied to a tree-level analy- 
sis of p(7, K + )A in section 4. We discuss the results of our work and highlight 
the important next steps in studying this problem. We present conclusions in 
section 5. 

2 Model Formalism 

We adopt a hadrodynamical framework for modelling kaon photoproduction 
on the nucleon. The model we use has been described in previous work [17], 
so we present only a brief summary here. The ^(7, K + )A reaction process is 
described by hadronic degrees of freedom using an effective Lagrangian. The 
contributioning diagrams are depicted in figure 1. Every intermediate particle 
in the reaction is treated as an effective field with associated mass, photo- 
coupling amplitudes and strong decay widths. Tree level Feynman diagrams 
contain the vector K*(892) and the axial- vector Ki(1270) t-channel mesons, 
as well as the usual Born terms. Two hyperon resonances, the Soi(1800) and 
P i(1810) are present in the ?/-channel. These ingredients constitute the so- 
called "background". 

In Ref. [18] we stressed the difficulties associated with parameterizing the 
background diagrams in ^(7, K + )A calculations and presented results for three 
plausible background schemes. Subsequent work [19] showed that the p(e, e'K + )A 
process is highly selective with respect to viable choices for dealing with the 
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Figure 1. Diagrams contributing to the p('y,K + )A process at the tree level. The 
upper row corresponds to the Born terms in which a proton is exchanged in the 
s-channel, a A or S + in the w-channel, and a K + in the ^-channel. The lower 
row shows the corresponding diagrams with the exchange of an excited particle or 
resonance. 

background diagrams. The model used here is the only one which we found to 
reproduce simultaneously the ^(7, K + )A and p(e, e'K + )A data. 

The finite extension of the meson-baryon vertices is implemented by the use 
of hadronic form factors [20,21,22]. We also note in passing that the validity of 
this approach is only likely to be reasonable for photon energies below about 
2 GeV, but its precise range of applicability remains to be established. 

For the purposes of the present work, we have defined several versions of 
our hadrodynamical model which correspond to different choices of s-channel 
resonances. For brevity, we subsequently refer to these as different model vari- 
ants, but the fact that they are based on the same hadrodynamical framework 
should be kept in mind. Each variant includes a "core" set of resonances, 
consisting of Sii(1650), Pn(1710) and Pi3(1720) which are well known to con- 
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tribute to the p(7, K + )A reaction [10,13,23]. The variant containing only these 
resonances is referred to as the "Core" variant. 

The other model variants each contain one other resonance of mass 1895 MeV. 
A D13 resonance was chosen to be compatible with our earlier work and be- 
cause it had been proposed in [10]. This in turn was motivated by the quark 
model of Capstick and Roberts [1] which showed that the Di 3 had the largest 
coupling strength to the lowest spin states. However, rather than relying on 
the predictions of a constituent quark model calculation, we have tried dif- 
ferent variants which contain a resonance of the same mass but alternative 
quantum numbers. If another resonance does contribute, the exact value of 
its mass is not likely to be important. Each of these variants is thus referred 
to by the character of the additional resonance: Su, Pn, P13 and D13 . Our 
previous analysis only used Core plus _D 13 resonances [14]. 

We therefore have set up our programme of calculations to address the ques- 
tions of whether an additional resonance is required in the description of the 
reaction (Core versus non-Core variants), and what character an additional 
resonance might have (Sn, Pn, P13 and D 13 variants). 

3 Analysis Procedure 

For each model variant, we performed a set of calculations to extract N* infor- 
mation from the reaction data. Each set consisted of 100 separate minimization 
calculations. The data sets employed in the fitting procedure included total 
cross-sections, differential photo-production cross-sections and recoil polariza- 
tions from the SAPHIR data set [11], and photon beam-polarization asym- 
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metries from SPring-8 [15]. In addition, we also used separated longitudinal 
and transverse electroproduction data from Jefferson Lab [16], as well as total 
electroproduction cross sections from a variety of older sources [24,25,26]. 

The coupling constants associated with the hadronic vertices in all the am- 
plitudes included in the model are free parameters in the fitting procedure. A 
description of their exact definition is given in [27], but we briefly list them 
here for completeness: coupling constants related to Born terms, gx+Ap and 
QK+Jip] vector and axial vector mesons in the t— channel, G V K *, G^*, G v Kl and 
G* K ; spin-| Y* resonances in the u— channel, Gy*Kp] spin-| N* resonances 
in the s— channel, Gn*ka] spin-| iV* resonances in the s— channel, G)j* KA , 
G\t*ka an d three off-shell parameters. Note that these constants are products 
of both a hadronic part and an electromagnetic part. Two hadronic cut-off 
parameters (one for all the Born terms, A& orn , and and for all resonances in 
the s— , t— and u— channels, A res ) are also free parameters. 

To compare the relative success of each model variant in describing the data, 
we require two separate but related procedures: a search for the best set of free 
parameters for each model variant, and a comparison amongst the different 
model variants. The first procedure is achieved by varying the free parameters 
and searching for a minimum x 2 statistic by fitting calculations to the data. 
This is a fairly standard procedure, but we remind the reader in passing the 
assumptions which are implicitly required: the data points being fitted must 
be independent (the value of one does not affect the others) , and error bars 
on the points represent the standard deviations of a Gaussian probability 
density. Both these requirements are likely to be approximately correct, and 
thus minimising ;\ 2 is equivalent to maximising a likelihood function. The 
likelihood represents the probability that the data would be measured, given 
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a particular set of free parameters. 

Furthermore, there is the assumption that no particular values of the free 
parameters are to be favoured prior to the fitting procedure. By imposing 
limits on the free parameters (as we do through physics constraints) this is 
not strictly true, but will be approximately true if the limits are sufficiently 
large. The whole procedure is then equivalent to maximising the probability 
that a particular set of free parameters is correct, given the experimental data, 
by varying the free parameters. 

Comparing different model variants needs to be done carefully, since each 
one may have a different number of free parameters. The approach we take 
implicitly employs "Occam's razor" by penalising model variants with more 
free parameters. 

3.1 Fitting 

Each model variant has at least 20 parameters which need to be extracted 
by fitting the calculations to the data. Traditional optimizing routines require 
a first guess at parameter values. Whilst some prior experience can be used 
to estimate the starting values, in a parameter space of this size it is very 
difficult to judge whether an optimum found by an optimizer is local or global. 
In previous work [17] a simulated annealing strategy was adopted, however we 
subsequently found [14] that using a genetic algorithm (GA) in combination 
with a traditional optimizer, MINUIT [28], offered many advantages. The GA 
is able to search a large region of parameter space and very quickly arrives 
at reasonable solutions, whereas a minimizer such as MINUIT, which uses 
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a Davidon-Fletcher-Powell algorithm, is able to take the solutions from the 
GA as starting points and find optima, provided the starting points are not 
too far from the optimum. This strategy thus plays to the strengths of both 
algorithms. 

Genetic algorithms are a class of search strategies known as evolutionary com- 
puting. A number of excellent texts on GAs exist (e.g. [29,30]), so we will only 
briefly sketch the strategy we have employed. The idea is to generate a num- 
ber of trial solutions randomly. In this implementation, each solution is an 
encoding of trial values of the free parameters as a string of real numbers 
{Aj}, where Aj represents the value of free parameter i. 

The collection of solutions is referred to as a "population". Each solution in 
the population is used to evaluate a function which determines its "fitness". In 
our case the fitness function is 

/{M= i + xW 

where x 2 {Aj} is the result of running the calculation of all the experimental 
observables and comparing with the available data. 

The population is then "evolved" in a manner analogous to biological evolu- 
tion. One or two solutions are selected from the current population, where 
solutions with greater fitness function values are selected preferentially, but 
not exclusively. When one solution is selected, it is subjected to a "mutation", 
where one or more of the free parameters are altered at random. 

When two solutions are selected, a new individual is created by "crossover" of 
the encoded parameters in each solution. In this implementation, crossover is 
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performed by a number of different functions, chosen at random, which fall 
broadly into two categories. The first is to swap one or more parameters from 
solution {Aj} to solution A, <-» /x, two obtain two "child" solutions. The 
second takes parameters from each solution and forms one child by assigning 
it new parameters {z/j} calculated by a weighted averaged: U{ = "^'t^' ■ 



A new fitness is then evaluated, and if it is better than the worst current 
fitness in the population, the new solution replaces the previous worst. This 
is often referred to as a steady state GA. The population therefore gradually 
migrates to one or more better points in parameter space, and if the routine 
is run for long enough, it will converge on one optimum. Whilst this optimum 
is not guaranteed to be global, GA research [30] has shown that "reasonable" 
solutions can be found very efficiently, even when no prior information has 
been used to constrain the free parameters. For instance, in our calculations, 
the best-of-generation solutions at the start of a run have x 2 values of order 
10 6 but this will be reduced to of order 10 after only 5000 evaluations of the 
objective function. If gradient methods were employed using such a random 
choice of initial starting values, it is highly likely they would be trapped in 
local minima. Note also that compared to the simulated annealing strategy 
used in previous work [17,18,19,31] we have reduced the number of function 
evaluations by a factor of 10. Whereas no one strategy can be optimum for 
all problems, the GA is likely to be highly efficient for problems such as the 
present one with many parameters and complicated \ 2 surfaces. The best-of- 
generation solutions at the end of each GA run can be used as starting points 
for a MINUIT minimization. 
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The complete strategy for analysis of each model was as follows: 

(1) 100 GA calculations were initiated. Each GA calculation used a popula- 
tion of 200 and was run for 5000 evaluations of the objective function. 
The limits for each parameter were deliberately chosen to define a large 
region of parameter space. 

(2) The best solution from each GA calculation was used as a starting point 
for a MINUIT minimization which was run for a further 5000 function 
evaluations. At this point very few, if any, of the MINUIT calculations 
had achieved convergence and the residual values of x 2 were too high; a 
more limited search space was required. 

(3) The parameters associated with each of the 100 MINUIT solutions were 
then examined. The set of solutions whose \ 2 values were within 1.0 of 
the current best x 2 were selected. A new, tighter set of limits for each 
parameter was defined from the range of parameter values exhibited by 
the chosen subset of solutions. This typically reduced parameter ranges 
by factors of 5-10. 

(4) 100 new GA calculations were initiated, the only difference from before 
being the smaller region of parameter space in which the search was 
performed. 

(5) As before, the best solution from each GA calculation was used to initiate 
MINUIT calculations. A few calculations had converged at this point. It 
was clear, however, that the calculations had not converged to the same 
optimum, so to investigate this, we proceeded to try to obtain convergence 
for all the calculations. 

(6) The MINUIT solutions were then re-inserted back into the GA popula- 
tions. Since the inserted solution was likely to be the "fittest" solution 
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in the population, this would cause the GA to concentrate more on a 
particular region of parameter space. Furthermore, the populations from 
each of the 100 GA calculations were mixed. This is a variation on what 
is known as an island population mixing scheme, and can be shown to be 
beneficial in cases where separate GAs are converging on different optima. 

(7) A third run of the 100 GA calculations was now performed. 

(8) A final MINUIT calculation for each of the 100 GA best solutions was 
done. At this point, all the calculations had converged. 

The result for each model variant was 100 solutions, each of which represented 
a converged minimization by MINUIT. As we noted in [14], we observed that 
each solution was unique, even although they had roughly similar x 2 values. 
In other words, there were at least 100 local minima in the x 2 surface. 

The fact that each calculation appears to find a different local minimum shows 
that the optimization problem may be ill-posed. The reason for this could be 
either that the reaction model lacks an ingredient which is necessary to ob- 
tain reproducible minimization results, or that the experimental data contain 
systematic uncertainties which render them inconsistent with each other, or 
that both the model and the data are problematic. From this observation, we 
caution that it is not possible to extract the magnitudes of coupling constants 
to high precision by fitting hadrodynamical models such as ours (and by ex- 
tension other similar models) to a relatively small number of photoproduction 
data points. 

What is observed however, is that the "best" of the calculations (i.e. those with 
the lowest x 2 values) tend to cluster around particular regions in parameter 
space. Our previous work [14] showed that this still resulted in considerable 
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ambiguities in the calculations of un-measured observables, and pointed to 
quantities such as photon asymmetry as being very useful in reducing the un- 
certainties in parameter values. The additional data used in the present work 
are the beam polarization data from the SPring-8/LEPS facility [15] and elec- 
troproduction data [16,24,25,26], and it appears they have affected a marked 
reduction in the regions of parameter space which give a good fit to the data 
(see section 4). This is a valuable observation. It shows that our original work, 
in agreement with others (e.g. [10]), was able to predict the measurements 
which would most significantly improve the comparison between data and 
theory. 



3.2 Model Variant Comparison 



Whilst the best fit for each model was obtained by minimising x 2 , this does 
not take into account the different number of free parameters in each model, 
or how small a region of parameter space corresponds to a reasonable fit. A 
full description of the framework we have employed is given in the appendix, 
but we sketch the salient points here. 

Comparing two model variants (A and B, say), we can evaluate the ratio of 
the likelihood functions: 

P(D\B)' K) 

where P (D\A) (P (D\B) ) is the probability that the data D would be ob- 
tained, assuming that variant A (B) were true. Under some reasonable assump- 
tions (as mentioned previously), minimising x 2 is equivalent to maximising the 
likelihood, but what is actually obtained could be written as P [D\Xf IP , A) 
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(where MP stands for "most probable"). This expression indicates that the 
maximum likelihood depends on the parameters of the model, Aj, whereas 
Eq.l requires there to be no dependence on the free parameters of either vari- 
ant. 

Using further approximations, we obtain for variant A the expression (with a 
similar one for variant B): 

yrN x\MP 

P (D\A) oc P (D\\ MP , A) „ ,7 ' , (2) 

where there are N free parameters, X™ ax and A™ n represent the maximum and 
minimum limits of search space for parameter A«, and 5Xf IP is the uncertainty 
in the most probable value of parameter Aj. The second factor in Eq.2 is 
often referred to as the "Occam factor", since it penalizes model variants with 
more free parameters. It is most easily interpreted as a product of N factors, 
each one a ratio of fitted error bar size to search space size for parameter Aj. 
Model variants whose parameter error bars are relatively large are thus more 
favoured by the data since it means that there is a larger region of parameter 
space which gives a good fit of the model to the data. The ratio in equation 
1 becomes not just a ratio of maximum likelihoods, but is multiplied by the 
ratio of Occam factors. 

Coupling constants associated with the Born amplitudes and the cutoff pa- 
rameters were free parameters in the fitting procedure, but tended to migrate 
to limits imposed by physical considerations such as SU(3) symmetry. How- 
ever, for all model variants these parameters were roughly similar, and so were 
"uninteresting" from a physical point of view. We therefore reduced the effec- 
tive number of free parameters after the fitting procedure by considering only 
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those parameters which related are to coupling constants of resonances in the 
s-channel. 

4 Results and discussion 

4-1 Fits to the available data 

The results from the fitting calculations are summarized in table 1. The factors 
quoted in the table are described as follows: 

(1) x 2 '- The value quoted is the lowest x 2 P er data point found for each model 
variant. 

(2) Number of free parameters: This is the number of s-channel reso- 
nance coupling constants which were varied in the fitting procedure, as 
explained in section 3. 

(3) Number of "best" calculations: For each model variant, a histogram 
of the number of calculations versus x 2 was produced. Fig. 2 shows three 
such histograms representing the results from the Pn, P13 and -D13 model 
variants. Remembering that each calculation is a converged MINUIT re- 
sult, it can be seen that there is a range of x 2 values for each variant. 
Examination of figure 2 reveals significant differences in the results ac- 
cording to model variant. For instance, Pn shows a large number of cal- 
culations at one value of x 2 ( a t the level of histogram binning), whereas 
P13 shows a much wider range in % 2 , but there is one calculation which 
has attained a much better % 2 than all the others, and _D 13 is somewhere 
in between the two. One can also see that Di 3 contains calculations with 
the lowest x 2 ■, but there are only a handful of them, whilst Pn has a large 
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Model Variant 


Core 


Sn 


Pu 


Pl3 


£>13 


x 2 


5.14 


4.47 


3.47 


3.75 


3.35 


Number of free parameters 


4 


5 


5 


6 


6 


Number of "good" calculations 


15 


1 


64 


1 


8 


Occam factor 


1.000 


3.278 


5.556 


0.018 


1.167 


Ratio of Posterior Probabilities 


1.000 


4.500 


12.571 


0.035 


2.786 



Table 1 

Results from the fitting calculations. The various factors are explained in the text, 
fraction of its calculations in the lowest x 2 "peak". The number of best 
calculations represents the number in the lowest x 2 "peak" for each model 
variant. 

(4) Occam factor: The Occam factors are also included to show how the 
various model variants are penalized for having either a larger number 
of free parameters, or which obtain good fits only in a smaller region 
of parameter space. They are normalized to the value obtained for the 
"Core" variant. 

(5) Ratio of Posterior Probabilities: The figures quoted are the result 
of applying the formulae derived in section 3 and the appendix. As with 
the Occam factor, they are also normalized to the value obtained for the 
"Core" variant. 

The "raw" x 2 results show that the P13 model variant gives the best fit to 
data (in agreement with other predictions [10]), but that the Pu and Pi 3 
variants also give comparable fits. However, once the Occam factor is taken 
into account, the ratio of posterior probabilities show that the Pu variant is 
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Figure 2. Histograms of the number of calculations resulting in x 2 values 

favoured by a factor of 4.5 over D\ 3 and a factor of 359 over Pi 3 . It is interesting 
to note that this mirrors the "intuitive" feeling we obtained when examining 
the distributions of \ 2 results shown in figure 2, i.e. the model variant with 
the greater number of low x 2 calculations (-Pn) is favoured over the one with 
very few (-P13). It is also interesting to note that S\\ is slightly favoured over 
-D 13 , even although the raw fit to data is not as good; this presumably reflects 
the fact that -D13 contains more parameters than Su. 

The main points we would like to conclude from this are: 



(1) The available data show evidence that, within our model framework, an 
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additional resonance is required. 
(2) The quantum numbers of the additional resonance are not possible to 
ascertain (given the amount of experimental data), but the present data 
set give greatest support to the hypothesis that it is a Pn. 

4-2 Comparison with the data used in the fitting procedure 

In order to show how well each model variant can fit the data, we present a 
selection of plots which show the comparison. For each variant, the parameters 
found by the fitting calculation which produced the best \ 2 value have been 
used. Our previous work [14] showed calculations using many different fitted 
parameter sets to make the point that a great deal of ambiguity remained 
after fits to data. However, now that we have included new data points in the 
fitting process, namely the photon beam asymmetry and electroproduction 
data, the ambiguity in parameters within the same model has been drastically 
reduced. We therefore restrict ourselves to plotting one calculation for each 
model variant, since the differences within a model variant are much smaller 
than the differences amongst the different variants. 

Figure 3 shows the comparison of all the model variants to the total ^(7, K + )A 
cross-section data, whilst figures 4 and 5 show typical differential cross-section 
and recoil polarization fits respectively. At this point we note that the four 
model variants which include an additional resonance are able to reproduce a 
"resonant bump" in the total cross-section, whilst the core variant cannot. 

The recoil polarizaton in figure 5 again shows that an additional resonance is 
required in order to reproduce the data especially at backward angles. Whilst 
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Figure 3. Total cross-section for the p(7, K + )A reaction. The data points are from 
[11]. The model variant calculations are: core set (solid line), Sn (dashed line), P\\ 
(small circles), P\3 (dot-dashed line) and D\z (solid plus circles). 



Figure 4. Fit of each model variant to differential cross-section data. The lines have 
the same meaning as in figure 3. 

a "bump" in the total cross-section was reproduced by Saghai [13] by altering 
an ingredient of the background processes, it is unlikely that the tuning of the 
background could result in a substantial non-zero recoil polarization as seen 
in the data. This observation may be a further indication that an additional 
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Figure 5. Fit of each model variant to the recoil polarisation data. The lines have 
the same meaning as in figure 3. 

resonance is required in this energy range. 

Figure 6 shows the calculations compared to the recent SPring-8 photon asym- 
metry [15]. It is clear that these data have greatly influenced the overall fit, 
since they tie down all the calculations at forward angles. Indeed, we predicted 
that this would be the case [14]. However, especially at the lower energy set- 
ting (1.55 GeV), it is still the case that there is a large difference in the model 
calculations over most of the angular range. Measurements over an extended 
range of angle would therefore be very useful to constrain further the fitting 
procedure. At 1.55 GeV, even an asymmetry averaged over the entire angu- 
lar range would discriminate between the Pn and Di 3 calculations, since the 
predictions are negative and positive respectively. 

Comparisons with electroproduction data are presented in figure 7. Longitu- 
dinal, transverse and combined differential cross-sections are shown. It is clear 
that agreement with the transverse cross section is the least good of all the 
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Figure 6. Fit of each model variant to photon beam asymmetry data. The data are 
from [15], and the lines have the same meaning as in figure 3. 

comparisons with data, and is probably responsible for the residual high val- 
ues of x 2 found in the fitting process. However, we note that there is virtually 
no difference amongst the model variants in the combined cross-sections, and 
only a small scale and form difference in the separated cross-sections. In other 
words, these electroproduction observables are relatively insensitive to details 
of the resonances in the reaction, and mirror our previous work [19] where we 
concluded that they are highly influenced by "background" processes. 

For completeness, and to summarize the results of our analysis, we show in 
table 2 the values of the resonance coupling constants which correspond to 
the best fit calculation for each model variant. In addition, we show graphi- 
cally the range and correlation of the coupling constant values in figure 8, for 
those coupling constants which are common to all models 1 . Error bars (where 

1 Note that in figure 4 of reference [14], there was an error in the normalisation of 
the Pi3(1720) and Z?i3(1895) coupling constants. The values used in the calculations 
of that article were correct. 
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Figure 7. Fits of model variants to electroproduction data. The separated longitudi- 
nal (L) and transverse (T) differential cross section filled square data are from [16] 
at W = 1.84GeV and 9* K = 0°, and the Q 2 = (photoproduction point) from [11] 
is shown as a comparison. The combined (T+L) differential cross-section data are 
from [24,25,26] at W = 2.15GeV and 6* K = 8°. The lines have the same meaning as 
in figure 3. 

shown) indicate the spread (standard deviation) of the values for the "best" 
calculations, defined in table 1, otherwise the spread is less than the size of 
the symbol. This plot indicates that the values extracted from the process are 
indeed highly correlated with the choice of additional resonance, as there are 
no discernable similarities in values and there is also a sizeable spread. 



4-3 Comparisons with other data 



A recent experiment using the CLAS detector at Jefferson Lab [32] has mea- 
sured the transferred polarization in the exclusive p(e, e'K + )A reaction. Figure 
9 shows the comparison of the model variants with this data set. We have not 
used this data in the fitting routine because the calculation of each point 
would require an averaging over the phase-space of each data bin, thereby sig- 
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Model Variant 


Coupling Constant 
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D 13 
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-8.1 x 10~ 2 ±3.6 x 1(T 5 


-4.9 x 1(T 2 


-7.2 x 10" 2 ±2.1 x 10~ 4 


-9.1 x 10~ 2 


-6.9 x 10" 2 ±2.2 x 10~ 4 


G[P U (1710)] 


-4.3 x 10" 2 ±9.4 x 1(T 5 


-1.1 x KT 1 


-4.4 x 10~ 2 ± 1.4 x 10" 3 


-1.1 x 10" 1 


-8.9 x 10" 2 ±3.5 x 1(T 3 


G( 1 )[P 13 (1720)] 


1.7 x 10" 2 ± 1.2 x 10" 5 


1.3 x 10~ 2 


2.2 x 10~ 2 ± 1.6 x 10~ 4 


2.2 x 10~ 2 


7.3 x 1(T 3 ±2.8 x 10" 4 


G( 2 )[P 13 (1720)] 


4.9 x 10~ 2 ±6.0 x 1(T 5 


-8.2 x 1(T 3 


3.6 x 10~ 2 ±6.9 x 10" 4 


2.4 x 1(T 2 


7.4 x 10~ 3 ± 1.8 x 1(T 3 


G[5n(1895)] 


- 
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-3.8 x 10" 1 ±5.5 x 10" 4 


- 
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G«[P 13 (1895)] 








2.1 x 10~ 2 
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-2.2 x 10 -1 ±5.7 x 10~ 4 
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-2.5 x 10~ 2 ±8.0 x 10~ 4 
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Figure 8. Depiction of the values of the main coupling constant found for each model 
variant. The meaning of the symbols are: circle - core set, square - Su, upright 
triangle - P\\, downwards triangle - P13 and star - -D13. 



nificantly increasing the computational time. Instead we show the predictions 
for each model variant, based on the extracted free parameters described in 
4.2. All the calculations show moderate agreement with the data, with the ex- 
ception of the P' x observable which was measured to be consistent with zero, 
where none of the model variants appear to be able to reproduce this be- 
haviour. It was suggested in [32] that this may be due to the production of ss 
quark pairs with anti-aligned spins. This is beyond the scope of the present 
hadrodynamical model. However, in common with the separated longitudinal 
and transverse differential cross-sections, there is relatively little sensitivity to 
details of which resonances are included in the reaction calculations, which 
indicates that this reaction is not best one to search for missing resonances. 
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Figure 9. Comparison with the double polarization electroproduction data of ref. 
[32]. The line conventions are the same as Fig. 3. 

4-4 Predictions for other observables 



We have already seen that a more comprehensive measurement of the pho- 
ton beam asymmetry would be highly sensitive to the choice of resonance in 
the reaction. With the prospect of double polarization photoproduction ex- 
periments [33], we present some predictions, based on the models described 
above and using the parameters found in the fitting procedures. The authors 
of [33] have proposed to measure beam-target, target-recoil and beam-recoil 
polarizations for both ^(7, K + )A and p(j, K + )Yj reactions at a photon beam 
energy of 1.5 GeV. 

Figures 10 and 11 show the results of the calculations using each of our models 
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Figure 10. Beam-recoil polarization prediction for linearly (top panels) and circularly 
(bottom panels) polarized photons. 

of beam-recoil and beam-target polarization for the p(-y, K + )A reaction only. 
Beam-target, target-recoil and beam-recoil polarizations are not independent, 
although in practice the measurement of all three would be desirable to limit 
systematic effects. All the observables show a high degree of sensitivity to the 
choice of resonance. One could therefore imagine that the measurement of a 
few selected points would be highly valuable in determining the dynamics of 
the reaction. 



5 Conclusion 

A framework for analysing p(7, K + )A reaction data sets has been presented. 
Using the foundation of a tree-level hadrodynamical model, we have shown 
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Figure 11. Beam-target polarization prediction for circularly (top panels) and lin- 
early (bottom panels) polarized photons. 

that the extraction of coupling constants by fitting to data is not trivial. The 
strong suggestion from our work is that fitting procedures should be carried 
out several times in order to establish how reliable the extracted numbers 
are. Carrying out the fitting procedure using a hybrid of a genetic algorithm 
and a traditional minimizer has been shown to play to the strengths of both 
strategies, and carrying out many calculations in parallel allows a thorough 
investigation of the character of the search space. Any conclusions from such 
procedures must take into account the uncertainty resulting from calculations 
which converge on different optima in parameter space. 

From the analysis shown here, we have reasonable evidence that an extra 
resonance is required to explain the ^(7, K + )A reaction data within our model. 
The study shows that the data currently favour the presence of an additional 
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Pn resonance of mass 1895 MeV, but that Su and Di 3 resonances are also 
weakly supported and cannot therefore be ruled out. We have not studied how 
combinations of additional resonances might improve the description. A firm 
conclusion awaits further data and a full-blown theoretical coupled-channel 
analysis to pin down the remaining uncertainty. There is a strong possibility 
that a GA may be of great value in the minimization required in a coupled- 
channel approach. 

We have shown that whilst the different model variants do a reasonable job of 
describing the p(e, e'K + )A data, they are relatively insensitive to the details 
of the contributing resonances. The polarization observables in photoproduc- 
tion reactions remain the best candidates for investigating the nature of any 
additional resonances. 

We are aware that further comprehensive data sets have been published [34,35]. 
We look forward to being able to use this to improve the quality of the infor- 
mation extracted in our approach. In addition, the inclusion of radiative kaon 
capture data [36] may provide further constraints on the values of resonance 
couplings. However, we note that in this reaction, which is crossed-symmetric 
to the p(7, K + )A studied above, coupled channel effects have been shown to 
be important [37,38]. 

This work was supported by the UK's Engineering and Physical Sciences Re- 
search Council, and the Fund for Scientific Research - Flanders. 
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Appendix: Model comparison 

We compare the relative success of the different model variants by evaluating 
the ratios of the maximum posterior probabilities [39]. To compare variants A 
and B the ratio, after application of Bayes Theorem, becomes 

P{A\D) _ P{D\A)P(A) 
P(B\D) ~ P(D\B)P(B) ' 

where P (A\D) is the maximum posterior probability for variant A, P (D\A) 
is the probability that the data would be obtained, assuming variant A to be 
true, P (A) is the a priori probability that variant A is correct and similarly 
for variant B. With no prior prejudice as to which variant is correct, we obtain 
the ratio of likelihoods: 

P(D\A) 
P(D\B)- 

The likelihood P (D\A) is an integral over the joint likelihood P (D, {A;} \ A), 
where {Aj} represents a set of free parameters: 

P(D\A) = fP(D,X i \A)d\ 

(3) 

= JP{D\\,A)P(X i \A) d\ 

where we have dropped the curly brackets for convenience. Model variants 
may have different numbers of parameters, and the multi-dimensional integral 
in equation 3 is the formal means of handling the problem. The function 
P (\i\A) is the prior probability that the parameters take on specific values. 
With no prior prejudice, we assume that each parameter Aj lies in the range 
ymm < ^ < X7 iax ? a nd we can write the prior as the reciprocal of the volume of 
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a hypercube in parameter search space (N-dimensional for N free parameters): 



P(Xi\A) 



nf i.K ax - K in ) 



To evaluate the volume of the hypercube in parameter space for each variant, 
the limits defined after the first GA plus MINUIT calculations were used, and 
in general they were different for each variant. 

The factor P (D\ Aj, A) in the integrand of equation 3 is not trivial to evaluate. 
If the errors in the fitted parameters, SXf IP corresponding to the maximum 
likelihood are multivariate Gaussian, it can be shown that 

P (D\\, A)^P (D\\f p , A) Det-i (J^j , (4) 

where H = — y 2 hi P (\\D, A) is the Hessian matrix, and is the inverse of the 
error matrix obtained in the fitting process which yields X MP . The explana- 
tions for this are given in standard texts on Bayesian data analysis [39]. 

Equation 4 can thus be written as 

P (D\A) oc P (D\\f /IP , A) — tt— (5) 

where the second factor is the Occam factor. We make a further simplification 
in our case. As mentioned in the main text, some of the free parameters, such as 
those respresenting coupling constants associated with Born amplitudes and 
cutoff parameters tended to migrate to limits. Apart from being physically 
"uninteresting", this also caused the determinant of the error matrix to be 
very sensitive to details of these fitted parameters, since the errors on each 
were very small. The physically interesting parameters are the ones which 
relate to resonances in the s-channel. 
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Rather than re-doing the whole minimization procedure with fixed Born ampli- 
tudes and cutoffs (which we deemed unnecessary due to the likely inaccuracy 
of the results from fitting to a limited data set), we took the errors on the most 
probable parameter values from the MINUIT calculations. MINUIT returns 
values which take into account the effect of correlations amongst parameters, 
so we applied the approximation 



The first factor in equation 5 can be evaluated up to a scale factor by 



Hence, by taking a ratio of the total likelihoods for two variants, all common 
factors will divide out, leaving a figure which represents the relative probability 
of two variants being correct. 
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